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canonical  momenta  invariants  perpendicular  to  the  beam  propagation  direction. 

The  linearized  Vlasov-Maxwell  equations  lead  to  an  exact  matrix  dispersion 

relation  which  is  valid  for  perturbations  about  general  beam  equilibrium 

Gq(pz)  and  includes  coupling  to  arbitrary  harmonic  number  (n)  of  the  fundamenta 

wiggler  wavenumber  kn.  No  a  priori  restriction  is  made  to  low  beam  density  (as 
2  2  2  u 

measured  by  u7c  L)  or  small  wiggler  amplitude  (as  measured  by  u>  /ckn  = 

/s  _  o  P  u  c  u 

eBAmc  kQ).  Moreover,  no  assumption  is  made  that  any  off-diagonal  elements 

in  the  matrix  dispersion  relation  are  negligibly  small.  A  detailed  numerical 

analysis  of  the  exact  dispersion  relation  is  presented  for  the  case  of  a  cold 

electron  beam  described  by  GQ(pz)  =  6(pz-pg).  It  is  shown  that  the  instability 

bandwidth  increases  rapidly  with  increasing  wiggler  amplitude  o^/ckg.  Moreover, 

except  for  very  modest  values  of  wiggler  amplitude,  it  is  shown  that  the  growth 

rate  calculated  from  an  approximate  version  of  the  dispersion  relation  can  be 

in  substantial  error  for  large  values  of  (k  +  nkg)/kg.  Preliminary  estimates  o 

the  influence  of  beam  thermal  effects  are  also  presented. 
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A  self-consistent  description  of  the  free  electron  laser  instability 
is  developed  for  a  relativistic  electron  beam  with  uniform  density  propa- 

*  /N  A  /S 

gating  through  a  helical  wiggler  field  Bu  =  -B  cos  kgzex  -  BsinkQzey.  The 
analysis  is  carried  out  for  the  class  of  solutions  to  the  Vlasov-Maxwell 
equations  described  by  f^fz.p.t)  *  ng6(Px)6(Py)G(z»pz,t)  where  Px  and  Py 
are  the  exact  canonical  momenta  invariants  perpendicular  to  the  beam  pro¬ 
pagation  direction.  The  linearized  Vlasov-Maxwell  equations  lead  to  an 
exact  matrix  dispersion  relation  which  is  valid  for  perturbations  about 
general  beam  equilibrium  Gg(pz)  and  includes  coupling  to  arbitrary  harmonic 
number  (n)  of  the  fundamental  wiggler  wavenumber  kg.  No  a  priori  restric¬ 
tion  is  made  to  low  beam  density  (as  measured  by  w2/c2k2)  or  small  wiggler 
amplitude  (as  measured  by  u>c/ckg  *  eB/ymc2kQ) .  Moreover,  no  assumption 
is  made  that  any  off-diagonal  elements  in  the  matrix  dispersion  relation  are 
negligibly  small.  A  detailed  numerical  analysis  of  the  exact  dispersion 
relation  is  presented  for  the  case  of  a  cold  electron  beam  described  by 
Go^Pz)  *  <$(pz  -  Po) .  It  is  shown  that  the  instability  bandwidth  increases 
rapidly  with  increasing  wiggler  amplitude  u)c/ckg.  Moreover,  except  for 
very  modest  values  of  wiggler  amplitude,  it  is  shown  that  the  growth  rate 
calculated  from  an  approximate  version  of  the  dispersion  relation  can  be 
in  substantial  error  for  large  values  of  (k  +  nkg)/kg.  Preliminary 
estimates  of  the  influence  of  beam  thermal  effects  are  also  presented. 
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1.  INTRODUCTION 


In  recent  years  there  have  been  several  theoretical  and  exper- 
imental  investigations  of  the  free  electron  laser  which  generates 
coherent  electromagnetic  radiation  using  an  intense  relativistic  electron 
beam  as  an  energy  source.  With  few  exceptions,  theoretical  studies  of  the 
free  electron  laser  instability  are  based  on  highly  simplified  models 
which  often  neglect  beam  kinetic  effects  and  coupling  to  higher  harmonics 

-V  ^  O 

of  the  fundamental  wiggler  wavenumber  k.Q,  or  make  use  of  very  idealized 
approximations  in  analyzing  the  matrix  dispersion  relation.  The  purpose 
of  the  present  paper  is  to  develop  a  fully  self-consistent  description  of 
the  free  electron  laser  instability  based  on  the  Vlasov-Maxwell  equations. 
The  final  matrix  dispersion  relation  [Eqi  (4&H  includes  all  beam  kinetic 
effects  and  coupling  to  arbitrary  harmonic  number  (n)  of  the  fundamental 
wiggler  wavenumber  (kri.  Moreover,  the  final  matrix  dispersion  relation 


makes  no  a  priori  assumption  that  any  off-diagonal  elements  are 
negligibly  small. 


The  present  analysis  assumes  a  relativistic  electron  beam  with  uniform 
cross-section  propagating  in  the  z-direction  through  a  helical  wiggler  field 
described  by  [Eq.  (2)] 

=  -B  cos  knze  -  B  sin  kn*e 
%  u  -vx  u  -vy 

A 

where  B  *  const,  is  the  field  amplitude  and  Aq  *  2ff/kQ  is  the  wiggler 
wavelength.  Moreover,  we  consider  the  class  of  exact  solutions  to  the 
Vlasov-Maxwell  equations  described  by  [Eq.  (12)]^ 


ffe(z,£,t)  -  n0<5(Px)<$(Py)G(z,pz,t) , 
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where  n_  =  const. ,  P  and  P  are  the  exact  canonical  momenta  invariants 
0  x  y 

[Eqs.  (6)  and  (7)],  and  spatial  variations  are  assumed  to  be  one-dimensional. 
A  detailed  analysis  of  the  linearized  Vlasov-Maxwell  equations  (Secs.  2.B 
and  3. A)  leads  to  the  matrix  dispersion  relation  [Eq.  (45)] 

U)2 

D>>  DLl<“>  '  1  + 


»  +  42) <«■>]) 


where  =  eB/ymc  is  the  relativistic  cyclotron  frequency,  co2  =  Arm^e  /ym 

-  2 

is  the  relativistic  plasma  frequency-squared,  yme  is  the  characteristic 
electron  energy,  is  a  constant  of  order  unity  [Eq.(43)],  the  suscepti¬ 
bilities  and  are  defined  in  Eqs.  (37)  and  (38),  D^(w)  is  the 

longitudinal  (electrostatic)  dielectric  function  defined  in  Eq.  (39),  and 
T 

Dn+^(oj)  are  the  transverse  (electromagnetic)  dielectric  functions  defined 
in  Eqs.  (40)  and  (41).  The  striking  feature  of  Eq.  (45)  is  that  the  dis¬ 
persion  relation  is  valid  for  arbitrary  harmonic  number  n.  (No  a  priori 
assumption  has  been  made  that  n  ■  +1).  Moreover,  the  effective  wavenumber 
variables  that  occur  in  the  various  factors  in  Eq.  (45)  are  k  +  nk^, 
k  +  (n+l)kp  and  k  +  (n-Dk^.  In  addition,  Eq.  (45)  describes  stability 
behavior  for  perturbations  about  general  beam  equilibrium  distribution 
Gq(pz),  and  no  a  priori  restriction  has  been  made  to  low  beam  density  (as 
measured  by  Wp/c  k^)  or  small  wiggler  amplitude  (as  measured  by  u^/ck^). 
Apart  from  the  linearization  assumption,  no  approximations  have  been  made 
in  deriving  Eq.  (45).  For  example,  we  have  not  assumed  a  priori  that 

DL(w)  'v  0  and  therefore  neglected  the  corresponding  term  on  the  right-hand 
n  — 


^ 1 


side  of  Eq.  (45).  The  latter  point  is  very  important.  In  Secs.  3.B  and 


3.C  we  present  a  detailed  numerical  analysis  of  the  exact  dispersion 

relation  [Eq.  (45)]  for  the  case  of  a  cold  electron  beam  described  by 

GQ(pz)  *  <$(PZ  -  Pq) •  The  exact  stability  results  are  then  compared  with 

the  approximate  results  obtained  from  Eq.  (45)  by  assuming  at  the  outset 
L  T 

that  Dn (uj)  ^  0  and  Dn_^(oj)  ^  0,  which  gives  the  reference  dispersion 
relation  [Eq.  (46)] 


dWu))d!  .(o>)  -  ~ 

n  n-i  l 


For  sufficiently  large  wiggler  amplitude  (as  measured  by  oj^/ck^) ,  it 
is  shown  in  Sec.  3.C  that  the  growth  rate  u).  »  Imto  obtained  from  the 
reference  dispersion  relation  [Eq.(46)]  can  be  in  substantial  error 


for  large  values  of  (k  +  nkg)/kg. 


For  completeness,  and  to  orient  the  reader,  it  is  useful  to  summarize 

here  the  interaction  wavenumbers  and  frequencies  pertinent  to  the  free 
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electron  laser  instability.  *  We  consider  a  cold  electron  beam  (Secs. 

3.B  and  3.C)  and  look  for  simultaneous  solutions  to  D^((i>)  *  0  and 

n 

T 

Dn_^(uj)  *  0.  Shown  in  Fig.  1  is  a  first-ouadrant  plot  (w  >  0  and 
k  +  nk^  >  0)  of 

U)  =  (k  +  nk0)yb  ±  “p/Y 

and 

a)  -  [c2(k  +  nkQ  -  kQ)2  +  w2  ]  X^2 


versus  k  +  nk^.  Note  that  there  are  generally  four  intersection  points. 
Solving  the  above  equations  simultaneously,  the  upshifted  wavenumbers 
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(intersections  with  k  +  nk^  >  kQ)  are  given  by 


<k+nk0>u  =  Y2kQ 


1+6, 


1  ^  !/^l 

1+2  -2-  JL 

~  ckn 


0  Wb  >  J 


+  yB  — E-k 

-  rPb  ckQ  ^0  ’ 


and  the  downshifted  wavenumbers  (intersections  with  k+nk^  <  k^)  are  given  by 


(k+hkQ^  =  y\ 


1-6, 


1  +  •)  — E_  _!_] 
ck.  r 


1/2 


l0  Y6b  J 


Cl) 


+  yg  — L.  ^ 

~  YPb  ckQ  0  * 


where  y  =  (1  -  $2)  and 

with  u>2/c2k2  «  1,  Y26u  ,  the 
p  U  b 

<k+-Vj  ■  <1+Bb>k2ko 


=  V^/c.  In  the  limit  of  low  beam  density 
intersection  wavenumbers  can  be  approximated  by 

■  V<1-eb>- 


and 


<k+”ko>d  ■  «-Vv2ko 


k0/(l+6b). 


For  sufficiently  large  y»  the  upshifted  wavenumber  (k+nk-)*  =  (1+3,  )Y2k„ 

u  u  b  U 

can  correspond  to  very  short  wavelengths  and  is  the  intersection  region 
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of  interest  for  the  free  electron  laser  instability.  *  Depending  on  the 
value  of  beam  density  (u2/c2kQ)  and  wiggler  amplitude  (c^/ck^),  however, 
it  is  important  to  note  that  the  dispersion  relation  (45)  may  also  support 
instability  in  the  long  wavelength  intersection  region  corresponding  to 
(k  +  nkg)*.  (See  Sec.  3.C  and  Fig.  3). 


/ 


2.  THEORETICAL  MODEL  AND  ASSUMPTIONS 


A.  General  Theoretical  Model 

The  present  analysis  assumes  a  relativistic  electron  beam  with 
uniform  cross  section  propagating  in  the  2-direction.  The  beam  density 
is  assumed  to  be  sufficiently  small  that  equilibrium  space  charge 
effects  are  negligibly  small  with 

E°=0  .  (1) 

'V, 

Moreover,  the  electron  beam  propagates  through  a  helical  wiggler 

12  5 

magnetic  field  described  by  *  ’ 

|0=-Bcosk0z|x-Bsink0zey  ,  (2) 

where  B=const.  is  the  field  amplitude,  X0*2n/kQ  is  the  wavelength, 
and  e  and  e  are  unit  Cartesian  vectors  in  the  plane  perpendicular 
to  the  propagation  direction.  The  vector  potential  associated 
with  Eq.  (2)  is  given  by 

^°-(S/k0)cosk02^x+(B/k0)sink0z^y  .  (3) 

It  is  also  assumed  that  the  beam  density  and  current  are  sufficiently 
small  that  the  equilibrium  self  magnetic  field  can  be  neglected  in 
comparison  with  ^ 

We  consider  perturbations  in  which  the  spatial  variations  are  one¬ 
dimensional  in  nature  with  3 /3  x-0-3 /3  y ,  and  3 /3  z  generally  non-zero. 
Introducing  the  perturbed  potentials, 

6$ (z, t) 

fiA^x.O-fiA^z.t^x^Ay^.t^y  » 


and 


the  electromagnetic  field  perturbations  6E  and  6B  can  be  expressed  in 
the  Coulomb  gauge  as 


6£(2*t)  =  -  C  k  6AX(Z>t)€x_  C  k  6Ay<Z*tJ«y  *  (4) 


6£=-  h  6K^^ix+  k  6Ax(*»t)i 


where  6£=-V64>-(l/c)  (3  /3  t)<5£  and  6£=Vx6£.  In  the  present  geometry, 

there  are  two  exact  single-particle  invariants  in  the  combined 

equilibrium  and  perturbed  field  configuration.  These  are  the  canonical 

momenta,  P  and  P  ,  transverse  to  the  beam  propagation  direction,  i.e., 
x  y 

P  =p  -  —  A^(z)-  —  6A  (z, t)=const. ,  (6) 

x*xcxcx 

P  =p  -  —  A^(z)-  —  6 A  (z,t)=const. ,  (7) 

y  y  c  y  cy 

where  A®(z)-(B/k0)coskQz  and  A°(z)=(B/k0)sink0z  [Eq.  (3)],  and  px 
and  p  are  the  transverse  mechanical  momenta. 

y 

The  potential  perturbations  6$(z,t),  <5A  (z,t)  and  6A  (z,t)  are 

x  y 

determined  self-consistently  from  the  Maxwell  equations 


\c  3 1  3  z  /  ; 

\c23t2  3Z2/  y  C  J 


d3P  Vx(fb-fb> 


d3p  Vy(fb-fb> 


— j  6$“4ire  f 
a  9  J 


d3P(fb-f°)  . 


where  f^(z,^)  is  the  equilibrium  (3/3t=0)  beam  distribution,  and  f^(z,jj, 
in  general  solves  the  nonlinear  Vlasov  equation 


^—  +  v  5 - e  |6E+ 

3 1  z  3z  'v 


vx(B  +6B)-\ 

%  %  ,v»  3 


In  Eqs.  (8)-(ll),  the  particle  velocity  ^  and  momentum  jg  are  related  by 


8 


l 

i 


2  22  1/2 

mrs/(1+s  ’/m  c  )  .  Moreover,  m  is  the  electron  rest  mass,  -e  is  the 

electron  charge,  and  c  is  the  speed  of  light  in  vacuo. 

For  present  purposes,  we  examine  the  class  of  exact  solutions  to 
Eq.  (11)  of  the  fomP 


fb(*,«.t)-n0«(Px)«(Py)G(z,pz,t)  ,  (12) 

where  n-=const. ,  and  P  and  P  are  the  exact  invariants  defined  in 
0  x  y 

Eqs.  (6)  and  (7).  Note  from  Eq.  (12)  that  the  effective  transverse 
motion  of  the  beam  electrons  is  "cold".  Substituting  Eq.  (12)  into 
Eq.  (11)  and  making  use  of  Eqs.  (6)  and  (7),  we  find  that  Eq.  (12) 
solves  Eq.  (11)  exactly  provided  G(z,pz,t)  evolves  according  to  the  one¬ 
dimensional  Vlasov  equation 

{fr+y2|7-fjii<2-t)lr}GU,p2,,;)‘0  •  (13) 

z 

where  H  is  defined  by 

H(z,t)**YTmc2-e6i)>(z,t)  .  (14) 


In  Eq.  (14) , 


2_  r  2  4  2  2  2/A0. x*  \2  2,A0JxA  .2,1/2 

YTmc  =  [m  c  +c  pz+e  (Ax+6Ax)  +e  (A  +6A  )  ] 


(15) 


is  the  particle  energy  for  P  =0=P  .  Moreover,  substituting  Eq.  (12) 

x  y 

into  Eqs.  (8)-(ll),  the  equations  describing  the  nonlinear  evolution  of 

the  potential  perturbations  can  be  expressed  as 

2 


lc23t2  3z2/  * 

I  2  „  2  .  2  i 

\  c  3 1  3z  / 


4irnQe 


me 


4irn„e 


6A  =  - 


me 


(A^+<5A  )  [ 
x  x  J 

^G-A° 

VT 

f^o)- 

(16) 

(A°+6A  ) 
y  y  . 

[  +*  G-A* 

1  yt  : 

11  Z  ‘of  • 

(17) 

- — -  6i|i»«-4TrenQ 
3z  * 


dPjIG-Co)  , 


(18) 
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where  YTmc  is  defined  in  Eq.  (15),  GQ(pz)  is  the  equilibrium  Q/3t  =  0) 
beam  distribution,  G(z,pz,t)  solves  the  nonlinear  Vlasov  equation  (13), 


2  ,  2  4^  2  2^  2-  2.1/2 

Ymc  =  (m  c  +c  p^+e  B/kg) 


is  the  energy  in  the  absence  of  perturbations  (<5A  =0=6A  ).  In  obtaining 

x  y 

Eq.  (19)  from  Eq.  (15),  use  has  been  made  of  (A^)2+(A^)2=B2/k2=const. 
fEq.  (3)]. 

Within  the  context  of  the  present  model,  Eqs.  (13)  and  (16)~(18) 
describe  the  exact  nonlinear  evolution  of  the  system  for  perturbations 
about  the  general  beam  equilibrium  Gq(p^).  In  this  regard,  we  note  from 
Eqs.  (13)  and  (14)  that  the  axial  force  on  an  electron  in  the  phase 
space  (z.p^)  is  given  by 

2 

F  =-  ~~  H - — x-  f-  [2 A°6A  +2A°6A  +(6A  )2+(6A  )2]+e  .  (20) 

z  3z  „  2  3z  x  x  yy  x  y  3z 

2y^nc  3 

That  is,  the  effective  ponderomotive  potential  is  proportional  to 

2A°6A  +2A°6A  +(6A  )2+(<5A  )2. 
x  x  y  y  x  y 


B.  Linearization  Approximation 


For  purposes  of  the  stability  analysis,  we  now  consider  the 
linearized  version  of  Eqs.  (13)  and  (16)-(18).  In  this  regard,  it  is 
useful  to  introduce  the  dimensionless  potentials  defined  by 

«♦  =  ’  (21) 
ymc 

~^~i  («Vi6Ay)  *  (22) 

ymc  3 

A°-  3^-2  (A°±iA°)-  exp  (±  ikQz)  .  (23) 

y  me  3  0 

_  2 

In  Eqs.  (21)-(23),  ymc  -const,  denotes  the  characteristic  mean  energy 
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of  the  electron  beam,  and 


w  =  - —  , 
c  ?mc 


is  the  relativistic  cyclotron  frequency  associated  with  the  wiggler 
field  amplitude  B.  Making  use  of  Eqs.  (21)-(24),  the  perturbed  force 
6F  [Eq.  (20)]  and  inverse  relativistic  mass  factor  [Eq.  (15)]  can  be 
approximated  by 

«FS~;»C2  {i  1  fj  [=xp(lk0z)SA.+exp(-ik0z)6i+]-  |j  4}  .  (25) 

and  „  . 

-2  ui 

1_  =  i  _  1  ^  [exp(ik0z)6A_+exp(-ikQz)6A+]  ,  (26) 

YT  Y  2  y  0 

2  2422  2*2  2  1/2 

where  -ymc  =(m  c  +c  p^+e  B  /k^)  [Eq.  (19)].  Substituting  Eqs.  (25) 
and  (26)  into  Eqs.  (13)  and  (16)-(18),  and  combining  Eqs.  (16)  and  (17) 
to  give  equations  for  6Ax±i6A^,  the  linearized  equations  for  6G(z,pz,t), 
6A+(z,t)  and  6<p(z,t)  can  be  expressed  as 

(ft  +  vz  h) 6G 

(27) 

=^mc2  {l  ^  y  h  ^exP (ik0z)<5A_+exp(-ik0z)6A+]-  ^  6*}  |^- 


dpz6C  , 


2  / 

/i  .  2  .  2  *  „  w  /  „  f  dp 

-  h)‘\~  (sv|  -f  G° 

cat  a  z  c  \  J  ' 

+  cf^  exp(ikQz)  |y/dPz  ~~  “  \  ^  [exp(ik0z)6A_+exp(-ik0z)6A+] 

-5s  co})'  <29> 


M  hjxtshat  rfiWi  i  t 


is  the  relativistic  plasma  frequency-squared,  and  64>,  6A+  and  u>c  are 
defined  in  Eqs.  (21),  (22),  and  (24).  In  the  limit  of  zero  wiggler 
amplitude  (<oc-K)) ,  we  note  that  Eqs.  (27)- (30)  give  the  usual 
uncoupled  electromagnetic  and  electrostatic  dispersion  relations  for 
perturbations  about  a  one-dimensional  relativistic  plasma  equilibrium. 
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3.  STABILITY  PROPERTIES 
A.  General  Dispersion  Relation 


For  the  equilibrium  configuration  considered  here,  the  electrons 

are  constrained  to  move  on  surfaces  with  P^=p  -(e/c)A^=0=p  -(e/c)A^=P^ 

x*x  x  y  yy 

2  24  2  2  *2  2  1/2 

and  ymc  =  (m  c  +c  p^+B  / Icq)  =const.  This  implies  that  the  axial 

momentum  p  is  a  constant  in  the  equilibrium  field  conf iguration. 
z 

Moreover,  the  corresponding  electron  trajectory  that  passes 
through  (z,pz)  at  time  t'=t  is  given  by 


z'=z+  —  (t'-t)  , 

ym 


(32) 


Pi=Pz. 


where  p  /ym«v  is  the  axial  velocity, 
z  z 

Without  loss  of  generality  we  expand  the  field  perturbations 
in  Eqs.  (27)- (30)  according  to 


^^“I^nexP  (k+nkg)z-iu>t]  , 
n 

exp(ikQz)6A_=^A^_1exp[i(k-»-nk0)z-iwt]  ,  (33) 

n 

exp(-ik0z)6^+-jA^+1exp[i(k+nk0)z-iwt]  , 
n 

where  XQ«2ir/k0  is  the  periodicity  length  of  the  wiggler  field.  Here 
we  assume  that  k  is  real  and  Imu  >  0,  corresponding  to  temporal  growth. 
A  completely  parallel  treatment  can  be  developed  for  spatially  growing 
perturbations  (Irak  <  0)  and  real  oscillation  frequency  u>.  Substituting 
Eq.  (33)  into  Eq.  (27)  and  integrating  along  characteristics,  we  obtain 
for  the  perturbed  distribution  function  6G(z,pz>t) 
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-  2  3C0 


<50ymc  - -  Ji(k+nkQ)exp(i(k+nkQ)z] 


z  n 


1 

2  ckQ  y  v“n-l ' "n+l/  Yn 


(A_  i+A_ 


(34) 


•L 


xj  dt ' exp [ i (k+nk  )v  (t '-t)-iwt ' ]  . 
U  Z 


Integrating  with  respect  to  t'  in  Eq.  (34)  with  Imw  >  0  gives 
_  2  3G0  (k+nkQ) 

6C=YmC  —  [  — ^  )y-  exp[i(k+nk0)z-iWt] 


O'  z 


(^n  2  ckQ  y  (An-l+An+l^)  ’ 


(35) 


2  24  2  2  ^2  2  1/2 

where  v  =p  /ym  and  ymc  = (m  c  +c  p  +B  /kn) 
z  z  z  U 

From  Eqs.  (28)- (30)  for  the  field  perturbations,  it  is  evident 


that  integrals  of  the  form  u^Jdp^SG  and  WpYjdpz6G/Y  are  required. 
Therefore,  comparing  with  Eq.  (35),  it  is  convenient  to  introduce  the 
effective  susceptibilities  x^(“>)  defined  by 


2-f 


(0)  -  2  2  f, 

X,j  (u))=Ymc  u)  Idp 


(k+nk0)3G0/3Pj 
Z  (t)~  (k+nkg)v 


(36) 


(1).  .  -  2  2-f 

Xjj  (td)«sYmc  UpYl 


dpz  (k-HnkQ)3G0/3pz 
Y  w-(k+nk0)vz 


(2)  ,  ,  -  2  2-2fdpz  0<+nk0)3G0/3pz 

*n  (u)=Y,nC  “pY  ~  'V(k4^kn'jv'“ 
*  y  U  z 


(37) 


(38) 


For  future  notational  convenience  we  also  introduce  the  longitudinal 
and  transverse  dielectric  functions  defined  by 


Dj;(w)=c2(k+nk0)2+x^O)(u>)  ,  (39) 

D^+1(u))=w2-c2[k+(n+l)k0J2-alW2  ,  (40) 

DT  .  (w)=u)2-c2  [k+(n-l)kf.]2-a1w2 
n-I  u  1  p 


(41) 
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where 


and 


dPz 

*H"t  w 

_-r  dP 

a3=^  J  ~r  Go(pz) 

J  y 


(42) 


(43) 


for  future  reference.  Note  that  the  constants  and  are  of  order 
unity  whenever  is  strongly  peaked  around  y=y. 

We  substitute  Eqs.  (33)  and  (35)  into  the  field  equations  (28)-(30) 
and  make  use  of  the  definitions  in  Eqs.  ( 36 )  —  (43) .  After  some  straight¬ 
forward  algebra,  we  obtain  the  matrix  equation  relating  An_i»  an<* 


-2 

0) 


-2 

u> 


_T  ,  1  c  ,  2.  (2).  1  c  ,  2,  (2), 

Dn+1+  2  2.  2  a3U)p+Xn  ^  '  2  2,  2  (a3“p+Xn  } 


c  k 


-2 


1  “c  ,  2,  (2). 

2  2.  2  (a3“p  Xn  J ' 
c 


c  kr 
G 

-2 

_T  .  1  “c  ,  2  (2). 

Vl+  2  “P  (°3VXn  >’ 

c  k0 


x(1) 

<*0  *n  • 


01 

c 

ck0  Xn 


(1) 


(44) 


The  condition  for  a  nontrivial  solution  to  Eq.  (44)  is  that  the  determinant 

of  the  matrix  vanish.  This  gives  the  dispersion  relation 

~2 


c  k0 


[Dn+l(uHDJ-l(a))] 


(45) 


xdx^  (o))]2-D^(o>)  ta3^+X^2)  («) ]  )  , 


which  determines  the  complex  eigenf requency  u>  in  terms  of  k+nkg,  k^, 

,  -2  ,  2.  2 
and  ui  /c  kA. 
c  u 

The  striking  feature  of  Eq.  (45)  is  that  it  is  valid  for  arbitrary 
harmonic  number  n.  (No  a  priori  assumption  has  been  made  that  n>*+l.) 
Moreover,  the  effective  wavenumber  variable  k'-k+nk^  occurs  in  every 
factor  in  Eq.  (45).  In  addition,  Eq.  (45)  describes  stability  behavior 


I 

r 

i' 

I 
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for  general  equilibrium  distribution  gq(pz)>  and  no  a  priori  restriction 

2  2  2 

has  been  made  to  low  beam  density  (as  measured  by  u^/c  k^)  or  small  wiggler 
amplitude  (as  measured  by  u^/ck^).  Since  Eq.  (45)  contains  no  approxima¬ 
tions  apart  from  the  linearization  approximation,  we  rexer  to  Eq.  (45) 
as  the  full  dispersion  relation  (FDR). 

In  circumstances  where  the  beam  density  is  low  and  the  wiggler  amplitude 

L  T 

is  very  small,  Eq.  (45)  supports  solutions  near  Dn(w)=0,  D^+^(w)=0 
T 

and  D  ,(<d)*0.  In  this  case.it  is  instructive  to  simplify  Eq.  (45) 
n-1 

7  T* 

near  the  simultaneous  zeroes  of  D  (oi)  =  0  and  D  ,  (w)=0.  (Here  we  assume 

n  n-1 

T 

Dn+^0.)  Equation  (45)  then  reduces  to  the  simplified  approximate  form 

-2 

DL(w)D*  .(*)=  [X(1)M12  •  (46) 

n  n— J.  z  z,  z  n 

c  kQ 

In  the  subsequent  analysis,  we  refer  to  Eq.  (46)  as  the  reference  dis¬ 
persion  relation  (RDR) . 

B.  Dispersion  Relation  for  a  Cold  Electron  Beam 

As  a  first  application,  we  consider  the  dispersion  relation  (45) 
for  the  case  of  a  cold  beam  equilibrium  described  by 


G0(Pz)“«(Pz"P0)  * 


(47) 


where  p^  is  related  to  the  mean  energy  ymc  by 


-  2 
ymc  ■ 


,  2  4X  2  2,  2  *2  ,.2.1/2 

(m  c  +c  pQ+e  B  /kQ) 


(48) 


It  it  straightforward  to  show  from  Eqs.  (47)  and  (48)  that  jdp2YnGQ(pz)“Yn 
so  that 


a^ayl  ,  (49) 

follows  from  Eqs.  (42)  and  (43).  Substituting  Eq.  (47)  into  Eqs.  (36)- 
(41),  the  effective  susceptibilities  and  dielectric  functions  for  a 
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cold  electron  bean  can  be  expressed  as 


x<0)M«-cW0)2 


2.-2 
u  h 
_ E _ 


[u>-(k+nkQ)Vb] 


2  ’ 


(1) . 2  {“Vc(k+t,k0>] 

Xn  (a,)-c(k+nk0)u,p  ———j  , 


X<2>  (w)  =  c(k+nkQ)(ij 


DL  (u>)  «  c2  (k+nk. ) 2 
n  u 


liu-(k+nk0)Vb] 

2  [2w6b-c(k+nk0) (1+B2) ) 
P  [u-(k+nk0)Vb]2 


1- 


2,-2 


[w-(k+nk0)Vb]‘ 


(50) 


(51) 


(52) 


(53) 


D^+1(w)=w2-c2 [k+(n+l)k0]2-w2  , 


(54) 


DT  .  (ai)=a)2-c2[k+(n-l)k_]2-a>2  . 
n-1  u  p 


(55) 


In  Eqs.  (50)- (55),  V^p^/yin  is  the  beam  velocity,  and  6b  is  defined  by 
Bb«Vb/c.  Substituting  Eqs.  (50)-(55)  into  Eq.  (45)  and  carefully 
ombining  and  simplifying  terms  on  the  right-hand  side,  the  full  dis- 
ersion  relation  (FDR)  for  the  free  electron  laser  instability  can  be 
t  ^pressed  as 


{ U-(k+nk0)Vb]2-Wp/Y2  }  {a,2-c2  ( (k+nk^-^J2-,,,2  } 

x  {a)2-c2  [  (k+nkQ)+k0 ]  2-cj2  } 

-2 

“c  2 r  2  2 .  ,2  2,  2  2, 

* - J~2  u  “c  <k+nk0>  -c  k0-lV 

r  k  **  ^ 

C 

X { [u- (k+nkQ) Vfa ] 2+2Vb (k+nkQ) [w- (k+nkQ) Vfe ] 
-c2(k+nk0)2/;2-Up  )  . 


(56) 


Note  the  exact  cancellations  that  have  occurred  on  the  right-hand  side  of 
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“"—2  2 

Eq.  (56).  Making  use  of  y  =1-8^,  the  full  dispersion  relation  (FDR) 
in  Eq.  (56)  can  be  expressed  in  the  equivalent  form 


{ [«- (k+nk0)Vb]2-w2/Y2 }  {a)2_c2  ( (k+nkjj)-^]2-^ } 

x  {u2-c2  [  (k+nk0)+k0]2-ui2 } 

. 2  (57) 

-  %>  -p[-2-c2(k4nk0)2-c2k2-w2] 
c  kQ 

x  [co2-c2  (k+nkQ)  2-to2  J  . 

Moreover,  substituting  Eqs.  (50),  (51),  (53),  and  (54)  into  Eq.  (46), 
the  reference  dispersion  relation  (RDR) ,  which  is  approximately  valid 
for  low  beam  density  and  very  low  wiggler  amplitude,  can  be  expressed  as 

{ [o»- (k-hik0) Vb ]  2-m2 /y2  }  (w2-c2  [  (k+nkQ) -kQ ] 2-o)2 } 

2  2  (58) 

1  “c  4  l“Vc(k-Hik0)]2 

“ 2  c2k2  u-a^w2  • 

2  2  —2 

Approximating  [u)-(k+nkQ)Vb]  «u)p/Y  in  the  denominator  on  the  right-hand 
side  of  Eq.  (57),  the  reference  dispersion  relation  (RDR)  becomes 


{[w-(k+nk0)Vb] 2-w2/y2 } {U2-c2 [ (k+nk0)-kQ] 2-m2 } 

.2 

“  I  ifl  V2[w8b-c(k+nk  )]2  . 
c  kQ 

C.  Numerical  Analysis  of  Dispersion  Relation 


(59) 


In  this  section,  we  summarize  the  results  of  numerical  studies 

of  the  reference  dispersion  relation  (59)  and  the  full  dispersion 

relation  (57)  for  a  broad  range  of  dimensionless  system  parameters 
2  2  2 

“c^ck0’  “p^c  kQ’  and  Y* 


Typical  results  are  summarized  in  Fig.  2  for  y«2  and  moderate  value 

of  wiggler  amplitude  (&c/ckg*0. 5)  and  beam  density  (ii)2/c2kg“0.4) .  For 

such  large  beam  density  and  wiggler  amplitude,  we  note  from  Fig.  2(a) 

that  the  growth  rate  u^-Imu)  obtained  from  the  reference  dispersion 

relation  (59)  can  be  in  substantial  error  for  large  values  of  (k+nkg)/kg 

This  follows  since  the  dashed  curve  in  Fig.  2(a)  [Eq.  (59)]  deviates 

significantly  from  the  solid  curve  (Eq.  (57)}  for  large  values  of 

(k+nkg)/kg.  On  the  other  hand,  from  Fig.  2(b),  both  the  reference 

dispersion  relation  (59)  and  full  dispersion  relation  (57)  give  very 

Similar  values  of  the  real  oscillation  frequency  ojr=Reu).  Moreover,  the 

oscillation  frequency  is  linearly  proportional  to  (k+nkg)  over  the 

entire  range  of  unstable  wavenumbers.  We  also  note  from  Fig.  2(a)  that 

the  maximum  growth  assumes  a  relatively  large  value,  [to. )  =0.3  ck_, 

i  max  0 

for  («c/ ckg,  uip/c2kg)«(0.5,  0.4). 

To  contrast  with  Fig.  2,  we  show  in  Fig.  3  a  plot  of  normalized 

growth  rate  n^/ckg  versus  (k+nkg) /kg  for  y=2,  moderate  beam  density 
2  2  2 

(o)p/c  kg-0.4),  and  small  wiggler  amplitude  (u>c/ckg=0. 1) .  In  this  case 

the  reference  dispersion  relation  (59)  and  full  dispersion  relation 

(57)  give  nearly  identical  growth  rate  curves.  Note  also  from  Fig.  3 

that  there  are  two  unstable  wavenumber  bands  [corresponding  to  two 

intersection  regions  of  the  u)  versus  k + nk„  curves  from  D^(w)  *  0  and 

U  n 

T  2  2  2 

D  ,(w)*0]  for  (&  /ckn,  w  /c  k_)«(0.1,  0.4).  Moreover,  the  maximum 
n-i  c  u  p  u 

growth  rate  assumes  a  relatively  small  value,  [u>. ]  =0.06  ck~,  for 

'  i  max  0 

the  small  value  of  wiggler  amplitude  chosen  in  Fig.  3. 

Shown  in  Fig.  4,  for  y"2  and  small  wiggler  amplitude  (&c/ckg-0. 1) , 

is  a  plot  of  normalized  growth  rate  u^/ckg  versus  (k+nkg)/kg  [Eq.  (57)] 

2  2  2 

for  several  values  of  normalized  beam  density  u^/ c  kg.  An  important 
feature  of  Fig.  4  is  the  fact  that  the  instability  bandwidth  remains 


2  2  2  2  2  2 

reasonably  narrow  for  u>  /c  krt  in  the  range  0.1  <  w  /c VZ  <  1.  This  is  in 

p  U  p  0  - 

contrast  to  Fig.  5  where  wi/ckQ  is  plotted  versus  (k+nk0)/kQ  [Eq.  (57)) 

2 

V 


—  2  2  2 

for  y=2,  moderate  wiggler  amplitude  (£  /ck  =0.5),  and  values  of  u>  /c  kt 

c  U  DU 


identical  to  Fig.  A.  The  rapid  increase  of  instability  bandwidth 

with  normalized  wiggler  amplitude  ai  /ck_  is  also  evident  from  Fie.  6 

c  u 

_  2  2  2 

where  the  growth  rate  is  plotted  for  y=2,  w  /c  k„=0.4  and  several 

P  0 

values  of  u>  /ck_  [Eq.  (57)]. 
c  u 

We  now  examine  the  case  where  y  is  large.  Figures  7-11  show 

plots  of  normalized  growth  rate  a^/ck^  versus  (k+nkg)/^  for 

Y*50  and  a  wide  range  of  values  of  normalized  wiggler  amplitude 

2  2  2 

u  /ck.  and  beam  density  w  /c  k-.  In  all  cases,  we  note  from 
c  u  p  0 

Figs.  7-11  that  the  instability  bandwidth  is  large.  This  is  in  contrast 

with  the  1ow~y  case  where  the  bandwidth  is  relatively  narrow  for  small 

values  of  S^/ck^  (Fig.  3).  We  also  find  (Figs.  7  and  8)  that  the 

growth  rates  determined  from  the  reference  dispersion  relation  (59) 

and  the  full  dispersion  relation  (57)  are  different  for  both  small 

wiggler  amplitude  (Fig.  8)  and  moderate  wiggler  amplitude  (Fig.  7). 

This  is  in  contrast  with  the  1ow-y  case  where  the  reference  dispersion 

relation  (59)  gives  a  good  estimate  of  the  growth  rate  for  small  values 

of  u^/ckg  (Fig.  3).  Comparing  Figs.  A  and  9,  for  small  wiggler 

amplitude  (ui^/ck^O.  1) ,  it  is  evident  that  the  growth  rate  is  considerably 

larger  for  large  values  of  y.  On  the  other  hand,  comparing  Figs.  5  and 

10,  for  moderate  wiggler  amplitude  (£oc/ckg=0.5) ,  the  growth  rate  is  only 

slightly  larger  for  large  values  of  y- 

2  2  2  -2 

Finally,  for  (k+nk_)  V  >>  w  /y  ,  we  expect  the  strongest  interaction 
*  U  D  p 

L  T 

between  the  longitudinal  and  transverse  modes,  D  (<u)=0  and  D  ,(w)=0, 

n  n-1 

to  occur  for  the  critical  value  of  wavenumber  given  in  Sec.  1,  i.e.,  for 


(  k+nlcQ }  u~  ( 1+®b  )  ^  2  k0 


(60) 
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For  6^=1  and  y=50,  Eq.  (60)  gives  (k+nkg)s=5000  kg.  For  y=50,  we  note 

from  Figs.  7-11  that  the  maximum  growth  rate  indeed  occurs  for  (k+nk_)/k..>5000. 

0 

D.  Influence  of  Beam  Thermal  Effects  on  Stability  Properties 


The  main  purpose  of  this  paper  is  to  develop  a  stability  formalism 
for  general  beam  equilibrium  distribution  gq(P2)  (Secs.  2.A-3.A)  and  to 
study  detailed  stability  properties  for  the  case  of  a  cold  electron  beam 
with  ^q^Pz)=<5 (Pz"Pq)  (Secs.  3.B  and  3.C).  For  completeness,  however, 
in  this  section  we  make  some  preliminary  estimates  of  the  influence  of 
beam  thermal  effects  on  stability  behavior.  As  a  simple  example,  we 
assume  a  beam  equilibrium  distribution  of  the  form 


Go(pz)=  T 


<pz"p0)2+A? 


(61) 


where  A^  is  the  characteristic  momentum  spread  about  the  mean  momentum  pQ. 
We  further  assume  that  the  momentum  spread  AQ  is  small  in  comparison 


with  the  directed  momentum  pg,  i.e.. 


A0  <<C  P0  * 


(62) 


Making  use  of  Eqs.  (36),  (39),  (61),  and  (62),  the  longitudinal  dielectric 
function  can  be  approximated  by 


D^(o,)-c2(k+nk0)2[l 


2,-2 

<yY 


[w-(k+nkQ)Vb+i  |k+nkQ  |  AQ/y3n]  ^ 


(63) 


-  7  7  A  7  7  7  ~7  .  7  \ I  7 

for  Aq  <<  Pq.  In  Eq.  (63),  yrac  =  (in  c  +c  p^+e  B  /kg)  is  the 


mean 


energy  [Eq.  (48)],  and  V  =p  /ym  is  the  directed  beam  velocity.  We 

b  U 

note  from  Eq.  (63)  that  the  fundamental  longitudinal  mode  obtained  from 
D^(w)=0  is  heavily  Landau  damped  by  thermal  effects  whenever 


(k+nk0)2  ^  . 

y  m 


(64) 
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It  is  useful  to  express  the  effective  momentum  spread  Ag  in 

,  .  _  .  .  _  2,24,2  2  2a2„  2,1/2 

terms  of  an  equivalent  energy  spread  Ay.  From  ymc  =(m  c  +c  p  +e  B  /kg)  , 

-  2  2  _ 

we  estimate  AY=PgAg/Ym  c  ,  or  equivalently  Ay/Y=8bAg/ymc.  Equation  (64) 
can  then  be  expressed  as 

c2(k+nk_)2(Ay/y)2>82  co2y2,  (65) 

U  bp 


for  Aq  «  Pq,  or  equivalently  for  Ay/y  <<  8^.  If  we  estimate  the  wave- 
number  of  interest  for  the  free  electron  laser  instability  by  (k+nkg)= 
(l+Bb)Y2k0  [Eq.  (60)],  then  Eq.  (65)  gives 


(Ay/y)  > 


«2  2 
Bu  i  w 


M*  1  m 

b  1  p 

2-2  22" 
(1+Bb)2  yVk2 


Equation  (66)  gives  an  estimate  of  the  energy  spread  required  for  heavy 

Landau  damping  of  the  fundamental  longitudinal  mode  obtained  from  D^(uj)=0. 

Whenever  Eq.  (66)  is  satisfied,  the  waves  are  heavily  damped  (linearly) 

and  the  free  electron  laser  instability  does  not  occur.  Equation  (66) 

constitutes  a  very  stringent  requirement  on  energy  spread.  For  example, 

2  2  2 

for  y=5,  6,  =1  and  u>  /c  kn*0.5,  Eq.  (66)  predicts  stabilization  for 
b  p  U 

(Ay/y)2  >  5><10  ,  i.e.,  for  a  fractional  energy  spread  Ay/y  of  about  7%. 

From  a  practical  point  of  view,  for  given  values  of  y  and  energy 

spread  Ay,  the  reverse  of  the  inequality  in  Eq.  (66)  can  be  used  to 

2  2  2 

determine  the  range  of  normalized  beam  density  w^/c  kg  required  for 
the  free  electron  laser  instability  to  occur,  i.e., 
w2  o  _2  (1+6.  ) 2 

■nft,  >  y  (Ay/y)  - 2 - -  •  (67> 

C  k0  6b 

Whenever  Eq.  (67)  is  satisfied,  it  is  important  to  note  that  the  inequality 
in  Eq.  (65)  will  provide  stabilization  of  the  instability  growth  rate 
for  sufficiently  large  values  of  wavenumber  (k+nkg) .  This  implies  a 
natural  tendency  of  beam  thermal  effects  to  limit  the  instability  bandwidth, 
which  should  be  contrasted  with  the  cold-beam  stability  results  in  Sec.  3.C. 
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4.  CONCLUSIONS 

In  this  paper,  we  have  developed  a  fully  self-consistent  description 
of  the  free  electron  laser  instability  based  on  the  Vlasov-Maxwell  equa¬ 
tions.  As  summarized  in  Sec.  2,  the  present  analysis  assumes  a  relativis¬ 
tic  electron  beam  with  uniform  cross  section  propagating  through  a  helical 

/V  A 

wiggler  field  B  =  -  B  cos  knz  e  -  B  sin  k  z  e  [Eq.(2)].  Moreover,  we 

u  ,\>x  u  r^y 

consider  the  class  of  exact  solutions  to  the  Vlasov-Maxwell  equations 

described  by  f,(z,p,t)  =  n_6(P  )6(P  )G(z,p  ,t)  [Eq.(12)).  A  detailed 
d  u  x  y  z 

analysis  of  the  linearized  Vlasov-Maxwell  equations  (Secs.  2.B  and  3. A) 
leads  to  the  exact  matrix  dispersion  relation  in  Eq.  (45).  The  striking 
feature  of  Eq.  (45)  is  that  the  dispersion  relation  is  valid  for  arbitrary 
harmonic  number  n.  Moreover,  Eq.  (45)  describes  stability  behavior  for 
perturbations  about  general  beam  equilibrium  distribution  Gq(Pz)»  and 
no  a  priori  restriction  has  been  made  to  low  beam  density  (as  measured 
by  oj^/c 2k^ )  or  small  wiggler  amplitude  (as  measured  by  u^/ck^).  In  Secs. 
3.B  and  3.C  we  present  a  detailed  numerical  analysis  of  the  full  disper¬ 
sion  relation  (Eq.(57)]  and  the  reference  dispersion  relation  [Eq.(59)] 
for  the  case  of  a  cold  electron  beam  described  by  gq(pz)  =  6(pz-Pq). 

Except  for  very  modest  values  of  wiggler  amplitude,  it  is  shown  in  Sec.  3.C 
that  the  growth  rate  =  Imto  obtained  from  the  reference  dispersion 

relation  can  be  in  substantial  error  for  large  values  o^  (k  +  nk^)/kQ. 

.// 

An  important  feature  of  the  stability  analysis  in  Sec.  f3,C  is  the  fact 

A 

that  the  instability  bandwidth  increases  rapidly  with  increasing  u^/ck^. 
Moreover,  the  instability  bandwidth  is  considerably  larger  for  large 
values  of  y  (Figs.  7  -  11)  than  for  moderate  values  of  y  (Figs.  2-6). 
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Finally,  in  Sec.  3.D,  we  have  made  preliminary  estimates  of  the  influence 
of  beam  thermal  effects  on  stability  behavior.  In  particular,  assuming 
that  Gq( pz)  is  given  by  the  Lorentzian  distribution  in  Eq.(61),  it  is 
found  that  waves  with  wavenumber  (k+nkQ)  =  (I+B^Y^q  are  heavily  Landau 
damped  whenever  the  fractional  energy  spread  Ay/y  exceeds  the  rather 
modest  value  in  Eq.  (66). 
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FIGURE  CAPTIONS 


L  T 

Fig.  1.  Intersection  of  dispersion  curves  (to)  =  0  and  ^  (w)  =  0. 

Fig.  2.  Plot  of  normalized  (a)  growth  rate  u^/k^c  and  (b)  real  oscillation 
frequency  ^/k^c  versus  normalized  wavenumber  (k  +  nk^/k^ 
for  y  =  2,  w^/ckg  =  0.5  and  io2/c2k2  =  0.4.  Solid  curves  are 
obtained  from  full  dispersion  relation  [Eq.  (57)1.  Dashed 
curves  are  obtained  from  reference  dispersion  relation  [Eq.  (59)]. 


Fig.  3.  Plot  of  normalized  growth  rate  UK/kQc  versus  (k+nkg)/kQ  for 

Y  =  2,  uWckg  =  0.1  and  U)2/c2k2  =  0.4  obtained  from  Eq.  (57) 
(solid  curves)  and  Eq.  (59)  (dashed  curves). 

Fig.  4.  Plot  of  normalized  growth  rate  tu^/k^c  versus  (k+nkQ)/k0  obtained 

A 

from  Eq.  (57)  for  y  *  2,  U^/ck^  =  0.1  and  several  values  of 

(u2/c2k2 
P  o 

Fig.  5.  Plot  of  normalized  growth  rate  u^/k^c  versus  (k+nk^/k^  obtained 

A 

from  Eq.  (57)  for  Y  *  2,  u^/ck^  =  0.5  and  several  values  of 

Fig.  6.  Plot  of  normalized  growth  rate  to^/k^c  versus  (k  +  nk^/k^  obtained 
from  Eq.  (57)  for  y  =  2,  u2/c2k2  =0.4  and  several  values  of 

“c/ck0* 

Fig.  7.  Plot  of  normalized  growth  rate  oj^/k^c  versus  (k+nk^/k^  for 

Y  •  50,  w  / ck0  =  0.5  and  u)2/c2k£  =  0.4  obtained  from  Eq.  (57) 

c  u  p  u 

(solid  curve)  and  Eq,  (59)  (dashed  curve). 


I 


Fig.  8.  Plot  of  normalized  growth  rate  (ib/k^c  versus  (k+nk^J/k^ 
for  Y  =  50,  L)c/ckg  =  0.1  and  U32/c2k.Q  =  0.4  obtained  from 
Eq.  (57)  (solid  curve)  and  Eq.  (59)  (dashed  curve). 

Fig.  9.  Plot  of  normalized  growth  rate  UK/kgC  versus  (k+nk^J/kg 

A 

obtained  from  Eq.  (57)  for  y  =  50,  ti^/ck^  =  0.1  and  several 

values  of  uj2/c2k2. 

P  0 

Fig.  10.  Plot  of  normalized  growth  rate  w^/k^c  versus  (k+nk^/kg 

A 

obtained  from  Eq.  (57)  for  y  -  50,  u^/ck^  =  0.5  and  several 

values  of  oj2/c2k2. 

P  o 

Fig.  11.  Plot  of  normalized  growth  rate  to^/k^c  versus  (k+nk^/k^ 

obtained  from  Eq.  (57)  for  y  =  50,  u2/c2k2  =  0.4  and  several 

A 

values  of  to  /ckn. 

c  u 
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